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Abstract 

We consider quantum correlations in a spin-1/2 open chain of N nodes with the XY Hamilto- 
nian using different bases for the density matrix representation and the initial state with a single 
polarized node. These bases of our choice are following: (i) the basis of eigenvectors of the fermion 
operators; this basis appears naturally through the Jordan- Wigner transformation (this represen- 
tation of the density matrix is referred to as the /3-representation), (ii) its Fourier representation 
(c-representation of the density matrix) and (hi) the basis of eigenvectors of the operators Ij z 
(the z-projection of the jth. spin, j = l,...,N). Although for the short chains (a few nodes) 
the qualitative behavior of the entanglement and the discord are very similar (the difference is 
quantitative), this is not valid for longer chains (N > 10). In this case, there are qualitative and 
quantitative distinctions between the entanglement and the discord in all three cases. We underline 
three most important features: (i) the quantum discord is static in the /3-representation, where 
the entanglement is identical to zero; (ii) in the c-representation, the concurrence may be non-zero 
only between the nearest neighbors (with a single exception), while the discord is nonzero between 
any two nodes; (iii) there is so-called "echo" in the evolution of the discord, which is not observed 
in the evolution of the concurrence. Using different bases, we may choose the preferable behavior 
of quantum correlations which allows a given quantum system to be more flexible in applications. 

PACS numbers: 
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I. INTRODUCTION 



The advantage of quantum correlations in comparison with classical ones rises a fast de- 
velopment of quantum information and communication devices. However, the problem of 
identifying of quantum correlations in a given system is not resolved yet. First, the so-called 
quantum entanglement l|-[5j was considered as a proper measure of quantum correlations. 
Later on, the quantum discord was accepted as a more adequate measure of quantum corre- 
lations p a|. In particular, the states with zero entanglement can reveal non-zero discord. 
Therefore it was noted Q that namely discord might be responsible for advantages of quan- 
turn information devices, in particular, for quantum speed-up. In general, the calculation of 
discord is very cumbersome optimization problem. In spite of intensive study of discord only 



very special cases have been treated analytically 
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4l3j. Nevertheless, namely these cases 



correspond to the reduced binary density matrix in spin- 1/2 chains governed by different 



15] or the initial state 



Hamiltonians with either the thermal equilibrium initial state 
with a single polarized node. 

A problem of a proper initial state is one of the fundamental problems in study of quntum 
correlations in different physical systems because of the technical difficulties of realization 
of a particular state. Most popular is so-called thermal equilibrium initial state which 
is most simple for realization. However, the initial state with a single exited node is more 

16 -El. 



relative, for instance, in quantum communication lines 
polarized node has been produced experimentally 22j. The evolution of quantum systems 
with this initial state at high temperatures has been studied, for instance, in 23), where the 
quantum echo has been found. 

Nevertheless, the problem of identification of quantum correlations is not resolved yet. In 
particular, the measure of quantum correlations depends on the basis which is taken for the 
density matrix representation. The reason is that, considering different bases, we involve 
different types of "virtual particles" . A possible way to avoid this ambiguity is suggested in 
ref. 24|, where the unitary invariant discord is introduced. This measure takes into account 
correlations among all possible "virtual particles". 

In this paper we consider the problem of preferable "virtual particles" (or the preferable 
basis of the matrix representation) from a different viewpoint. Instead of collecting the 
correlations among all "virtual particles" (like in the unitary invariant discord 24J) we con- 
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sider separately the quantum correlations among three types of particles, namely, among the 
fermions, which appear in a spin- 1/2 system with the nearest neighbor interaction under the 
Jourdan-Wigner transformation 25| (we call them as the ft- fermions) , among the fermions 
which are the Fourier representations of the /3-fermions (the c- fermions), and among the 
spin-1/2 particles in the basis of eigenvectors of the operators Ij Z (z-projection of the jth 
spin, j = 1, . . . , N). Each type of particles corresponds to the proper basis of the density 
matrix representation. We show, that the choice of the basis is a significant (may be even 
dominant) factor in identification of quantum correlations. 

This paper is organized as follows. The evolution of the spin-1/2 open chain of N nodes 
under the XY Hamiltonian with the nearest neighbor interactions is derived in Sec JTTI using 
the Jourdan-Wigner transformation. The general formulas for the discord and the concur- 
rence are discussed in Sec JIIII The comparison of the discord and the concurrence in different 
bases of the density matrix representation is given in Sec JIVI The summary of our results is 
represented in SecJVl 



II. SPIN-1/2 CHAINS WITH A SINGLE INITIALLY POLARIZED NODE 

In this paper we study the quantum correlations in the one-dimensional open spin-1/2 
chain of N nodes governed by the XY Hamiltonian in the approximation of nearest neighbor 
interactions, 

N N-l 

H = 0J 2^ hz + D / ](IjxI(i+l)x + Iiyl(i+l)y), (1) 
i=l i=l 

where Uq is the Larmour frequency in the external magnetic field, D is the spin-spin coupling 
constant between the nearest neighbors and Ij a (i = 1, . . . , N, a = x,y,z) is the ith spin 
projection on the a-axis. We chose the initial state of this chain with the single polarized 
jth node (1 < j < N ) at arbitrary temperature, i.e. 

Po = — = w [1 + 2I jt tanh |J , Z = Tr(e^) = 2 N cosh |, (2) 

where (3 = -^r, h is the Plank constant, k is the Boltzmann constant, and T is the tem- 
perature of the system. The evolution of the density matrix is described by the Liouville 
equation -£ = —i[H,p] for the density matrix p. The solution to this equation reads: 

p{t) = e- UH Po e ltH = ^e~ lHt (l + 2I ]Z tanh ^)e mt . (3) 

3 



To solve the Liouville equation, we diagonalize the Hamiltonian ([I]) using the Jordan- Wigner 



transformation method 



25] 



H = £ k^f3 k - ^Nco , e k = D cos(k) + 



wo, 



(4) 



where the fermion operators ftj are defined in terms of other fermion operators Cj by means 
of the Fourier transformation 



A' 



0k = ^<7fc0')« 

3=1 

and the fermion operators Cj are defined as j^] 



(5) 



Cj = (— 2) j I\ z l2 Z ■ ■ ■ I z {j-i)Ij 



Here 



9k{j) 



1/2 



sin(fcj), k 



n = l,2,...,N. 



(6) 



(7) 



N + 1J y Jil N + 1 

It may be readily shown that the projection operators Ij Z can be expressed in terms of the 
fermion operators Cj as 



/• =c + c- 



(8) 



Then the density matrix ([3]) can be transformed to the following form 



23] 



Pit) 



1 — tanh % tanh 
+ 



2 N 



2 n- 



-it[e k -e k ,) 



k,) 9kU)9k'U)PtPk'. 



(9) 



kM 



where we use the identity 



(10) 



It is interesting to note that the quantity 



for the pth node does not depend on /3 and 



coincides with the result of ref. [23j obtained for the high temperatures only ((a) = Tr{pa} 
for any operator a): 

2 



4 



(ipz)(t) 

(i jz )(o) (N + iy 



e sin(fcj) sin(fcp) 



(11) 



We will study the quantum discord Q nm and the concurrence C nm (as a measure of en- 
tanglement) between any two particles. Since Q nm = Q mn and C nm = C mn , hereafter we 
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take m > n without loss of generality, n = 1, . . . , N. As was mentioned in the Introduc- 
tion, we consider the open spin- 1/2 chain using its density matrix representations in three 
different bases: (i) the basis of eigenvectors of the fermion operators (3j (we refer to this 
representation of the density matrix as the /3-representation), (ii) the basis of eigenvectors 
of the Fourier transformed fermion operators Cj (the c-representation of the density matrix), 
and (iii) the basis of eigenvectors of the operators Ij Z (j = 1, . . . , N). Of course, consid- 
ering the correlations between the nth and the mth nodes in all these cases we consider 
the quantum correlations between physically different particles. Thus, using either the (3- 
or c-representation we consider the correlations either between the nth and mth /3-fermions 
or between the nth and mth c-fermions, while using the basis of eigenvectors of the opera- 
tors Ijz, we elaborate quantum correlations between the natural spin-1/2 particles. This is 
the reason of qualitative and quantitative difference among the discords and entanglements 
found in all these cases. 

First step in calculation of either discord or entanglement between the nth and mth nodes 
is the construction of the reduced density matrix with respect to all nodes except for the 
nth and mth ones. Below, we calculate the reduced density matrices in three different bases 
mentioned above. In all cases we use notations 

|nm) = {|00),|01),|10),|ll)} (12) 

for the basis vectors, where n in the vector \n) means the different filling numbers in the 
fermion bases or the excited spin (n = 1) and the ground state spin (n = 0) in the basis of 
eigenvectors of operators Ij Z . 



A. The reduced density matrix in the basis of eigenvectors of the fermion oper- 
ators (3j 

We use the superscript to specify the (3- representations of quantum operators. Let 
us reduce density matrix with respect to all fermions except for the nth and mth ones 
obtaining: 

1 tnnti tnnti ^ 

pL = \~ =±(£U) + 9lU)) + =± £ e-^-^ 9k (j) 9k/ (m^- (13) 

k,k'=n,m 
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Using the basis (1T2]) we find 



where 



rnm 



( f + //3 1 Tfi 







f + 7/3 7/9 
^00 ~ u mm °mn 



'J 



nm 





\ 





f + P n 



jS, / 



(14) 



J; 



1 tanh 2 



00 



1 I „2 



(^(j')+^(j'))> 



(15) 



It is obvious that 



tanh 2 



p _ tanh 2 2 



(16) 



which does not depend on the time t. 



B. The reduced density matrix in the basis of eigenvectors of the fermion opera- 
tors Cj 

In this case we write the density matrix (J2j) as 

1 farih S. tnnh ^ 

P {t) = 2N 2 + ==f E ^-^o^^awco^- (17) 

We use the superscript " c" to specify the c-representations of quantum operators. Reducing 
the density matrix p with respect to all fermions except for the nth and mth ones, we obtain 

1 tanh — 

P c nm = - A - E E e-«"'-*'' ) 9kU)9vV)9k{t)gv{t) + (18) 

fc,fe' l=n,m 

tanh I 



E E e^^^^o^jkW^co^'- 



2 

£:,£;' l,l'=n,m 
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The matrix form of the reduced density matrix p c nm in the basis (JT2]) coincides with eq. 
up to the replacements J% m — > J^ m with 

1 tanh 7r 
4 4 



7 C 



\ \ " e -i*(efe-£fc') 
fc,fc' l=n,m 



C ) 
mm) i 



T° - 

nm 

It is obvious that 



- --( r + /■ 



fc,/c' 



k '9k(j)9k'(j)9k(n)g k '(m). 



tanh | 



tanh£ 



(19) 



(20) 



where 5j n is the Kronecker symbol. Eq.f l20|) means, in particular, that J^ n < 1 



C. The reduced density matrix in the basis of eigenvalues of operators Ij z , j = 
l,...,N 

In this case we write the density matrix as 

Pit) = 2 7 2 E e-^-^ 9k (j)9Aj)9 k (l)9An x (21) 

{ — 2) l + l 2 I lz . . . I(l_l) z I^Il z . . . J(p_l) 2 /;7. 

We use the superscript "spin" to specify the matrix representations of quantum operators 
in the basis of eigenvectors of the operators Ij Z . Reducing the density matrix p with respect 
to all nodes except for the nth and mth ones, we obtain that the reduced density matrix is 
non- diagonal only if m = n ± 1: 

i tanh — 

aCi> = 7 - —r± E E *-** h - 9y) 9k{i)gw{j)9k{l)gv{t) + (22) 



4 4 

fc,fc' i=n,?i+l 



tanh | 



fc,fc 

9k{n)g k >{n)I+I~ + # fe (n + l)# fc /(n + 
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Otherwise (m/n±l) the density matrix is diagonal: 

P s r = ^ - — ^ £ E ^^g h {j)9M9k(t)9^l) + (23) 

l=n,m 

__JL £ e-*^-^) 5fc (i)^(i) [9 k (n)g kl (n)I+I- + g k (m) 9k ,(m)I+I- 

fc.fc' 

so that the discord and entanglement are both zero. The matrix form of the reduced density 
matrix p s ^+i) 1X1 ^ e basis f|T2|) coincides with eq.f lT4|) up to the replacements J% m — > J^ 1 
with 

jspin jc j spin jc jspin jc jspin jc (OA\ 

J 00 — J 00' J nn ~ J nm J ( n+ l)( n+ l) — ^( n +l)( n +l), J n {n+l) ~ J n{n+1)- 



III. GENERAL FORMULAS FOR DISCORD AND CONCURRENCE 

Since the reduced bi-particle density matrices in all three cases have the same matrix 
form (Tbf|) . we derive the formulas for both the concurrence (as a measure of entanglement) 
and the discord for the density matrix ffT4|) omitting the superscript (3. We will use the 
relation 

I Jmn | JnnJmm (25) 

which is valid in all three cases considered in Secs lII ATlII CI 
A. Concurrence 

We characterize the ent ang lement by the Wootters criterion in terms of the concurrence 



26|,|27j). According to 



26 



271 ]. one needs to construct the spin- flip density matrix 



P(nm)(r) = (a y ® ff,)(pW)'( T )( a( ® a y ), (26) 



where the asterisk denotes the complex conjugation in the basis (|12[) and the Pauli matrix 
cr y = 2I y . The concurrence for the density matrix p( nm )(r) is equal to 

C = max(0, 2A — Ai — A 2 — A 3 — A 4 ), A = max(Ai, A 2 , A 3 , A 4 ), (27) 
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where Ai, A2, A3 and A4 are the square roots of the eigenvalues of the matrix product 
P{nm){ T )P{nm){ T ) ■ For the density matrix p given by eq. ffT4"]) we have 



A = A 



1 . 4(t/ mm J nn ) -|- yj ' JmmJnni (28) 



A2 . 4(i/ mm J nn ) ~\J ' JmmJnm 



A3 — A4 — — y 1 — 4(J mm + J nn ) 2 . 



Substituting eqs. (l28|) into eq. (!27I) we obtain 

Cnm max ^0, 2 \J J mm J nn 2 4(t/ mjn -|- Jnnj^j . 



(29) 



B. Discord 



As far as matrix (j!4p is a particular case of the X-matrix, we will use the results of ref. 



lgj , where the discord for the X-matrix has been studied. We remind the basic formulas 
for calculation of the discord Q between two particles, which are called the subsystems 
A and B 12J. Recall, that the superscript B means that the projective measurements are 
performed over the subsystem B. The discord is introduced as a difference between the total 
mutual information J(p) encoded into the system AB and its classical part C B , 

Q^ = I{p)~C^ B \ (30) 

The total mutual information I reads 

Hp) = S(pP>) + S(p^) + J2 Ai log 2 Aj. (31) 

i=o 

Here p^ and p 1 - 5 - 1 of the reduced density matrices of the subsystems A and B, S(p^) and 
S(p^) are the von Neumann entropies of the subsystems A and B, 

S(P {A) ) = -\ (\og 2 {\ - JlJ + 2J mm log 2 , (32) 



1 A ,1 .9 x „ T , 1 + 2 J, 



S(P (S) ) = -o l0 g 2 ( 7 - + 2Jnnl0g 2 



■fin 



nn 



Xj are the eigenvalues of the density matrix f lT4"j) , 

Aq,1 = t(1 — 2(J mm + J nn )), A2,3 = ^(1 + 2(J mm + J n n)) ■ (33) 
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As far as all eigenvalues of any density matrix must be non-negative, eqs. (!33|) mean 

Jmm Jnn — 1/2- (34) 



The classical counterpart reads 12 



= S(pM) - min (poSo + PiSi)- (35) 

»7={0,1} 

Here the conditional entropies Si (i = 0, 1) are 

1 - 0W , 1 - 0« 1 + 0(0 l + 0(0 
Si = 2 — g2 — 2 2 — 2 — 2 — ' ^ 



0(0) 



2\/ Jmm(Jmin 1) Jnn) 2\J J mm {^J mm (jP 1 ) J nn ) 



1 + 2?7 J nn 1 - 2r?J nn 

and the populations p iy i = 0, 1 are 

Po = ^(l + 2r/J nn ), pi = -(1 -2r)J nn ), (37) 



where we introduce the parameter rj instead of k used in ref. [l2|, k = ^y 2 . It might be 
readily demonstrated that the minimum in eq.( |35l) corresponds to i] = 0. For this purpose 
we show that the derivative of the function 

f(v)=PoS +PiS 1 (38) 

with respect to the parameter 77 is positive over the interval < 77 < 1. In fact, this derivative 
reads after some transformations: 

™ = < a + »^ \l%t-% + (39) 

, a \ -i l + 2r]J nn + 2d l + 2r]J nn 

{(3 - a) log 2 — — — + 2a log 2 



1 - 2r]J nn - 2d 1 - 2r]J nn 

"nn a V^mmJun , f~r 7~t 7 2 -r \ t \" 



Collecting terms with a and /3 we obtain: 



f'M - a log ^ 2 (lWnn) 2 -(l-4,Vl) 2 (1 + 2^-4^ 

/ W - ^ 4^(1 - 2* 7 J TM ) 2 - (1 - VJD* + ^° g2 (1 - 2d? - 4*?JL' 1 } 

Both terms in the RHS of eq. ()4Up are nonnegative. Thus we conclude that f'(f]) > 0, i.e. 
f(j]) is increasing function over the interval < 77 < 1. Consequently, the function f(rj) 
takes the minimal value at the boundary point rj = 0. 
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As a result, the expression for the classical correlations (I35I) reads: 

C l j = ^log 2 _ + J mm log 2 + (41) 

mm ' ^"Jmm 

f~f j t T , 1 "I - 2 \j JmrrSyJmm Jnn) 

V Jmm\Jmm "T ^nn) !Og 2 ~ ~ ; , . 

1 2 •y/ Jmm \ Jmm -Jnn ) 

Now the discord may be calculated by formula ( 13T)j) using eqs. ()31H35j) and fl4"Tj) : 

= -^((l-2J nn )log 2 (l-2J nn ) + (l + 2J nn )log 2 (l + 2J nn )- (42) 
(1 2J mm 2J nn ) log 2 (l 2J mm 2J nn ) 
(1 "I - 2t/ mm -|- 2</ nn ) log 2 (l -|- 2J mm -\- 2J nn ) -\- 

(1 2 \/ Jmm^yJmm ~t~ Jnn)) log 2 (1 2 ^/ ' Jmm(Jmm "I" Jnn)) 
(1 "I - 2^Jmm(^mm "I - Jnn)) log 2 (l "I - 2^/ ' Jmm^Jmm Jnn)) ) • 



Note that the discord (obtained by means of the Neumann type measurements over 
the subsystem A) differs from the discord Q^ 5 - 1 by the replacement n -H- m. As far as 
7^ m general , we define the discord Q nm as follows [28j | : 



g Bm = min(g<£,ga)- ( 43 ) 

IV. DISCORD AND CONCURRENCE IN DIFFERENT DENSITY MATRIX 
REPRESENTATIONS OF THE OPEN SPIN CHAIN WITH N > 10 NODES 

The problem of description of those quantum correlations which are responsible for advan- 
tages of quantum computations is not resolved yet. Therefore it is reasonable to characterize 
different types of quantum correlations in a quantum system, i.e. not only those which cor- 
respond to the real particles, but also correlations between the "virtual particles". This was 



the basic motivation for introduction of the unitary invariant discord in ref. 24] as a way 
to count all quantum correlations between all possible "virtual particles". The geometric 
measure Q G of the unitary invariant discord is representable in terms of eigenvalues \ of 
the considered density matrix, 



Q G = ^ , in) 



so that it achieves the maximal value 1 for a pure state and is zero for the matrix with all 
equal eigenvalues. In the case of bi-particle density matrix f lT4|) . N = 2, with eigenvalues 

11 



(|33|) we obtain 





(45) 



The significant value of Q suggests us to look for such basis where the quantum correlations 
are valuable and/or the evolution of quantum correlations is suitable for applications. 

In the present paper we investigate the quantum correlations in a spin- 1/2 chain of N 
nodes using three different matrix representations, i.e. we consider the same chain using (i) 
the /3-representation introduced in Sec JII A\ (ii) the c-representation introduced in Sec lII Bl 
and (iii) the representation in the basis of eigenvectors of the operators Ij Z , j = 1, . . . , N, 
discussed in Sec iH CI Below we study the quantum correlations in all three bases and demon- 
strate the significant qualitative and quantitative distinctions among the discord and the 
concurrence in all three cases. It will be illustrated that the discord in the /3- and c- 
representations of the density matrix reveals the new properties of the considered quantum 
system in comparison with the usual discord in the basis of eigenvectors of the operators 



A. Discord and concurrence in a system of /3-fermions 

First, we consider the (3- representation of the density matrix (see Sec JII Ap and calculate 
the discord and the concurrence between the nth and mth fermions. Formula ( 12~T)) for the 
concurrence reads in this case: 



Simple estimation at (3 — > oo shows that the concurrence may be positive only for N < 4, 
which is not interesting for us. In the case of longer chain, C^ m = for any nodes n and m. 

The behavior of the discord is completely different. Using formulas ( I42p and ( I15p we 
obtain that the discord is non-zero in general and does not evolve in time. The discord Q n ,n+i 
between the nearest neighbors versus (3 = for N = 21 is depicted in FigJTJ Emphasize 
that the discord is non-zero not only between the nearest neighbors, but between the remote 
nodes as well, which is not valid in general, see, for instance, Sec JIV Cl The dependence of 
the discord Q^ n+l , I = 1,2,3,4, on n for iV = 21 and large values of (3 (f3 = 10) is shown 
in Figj2]for different initially polarized nodes j = 1,6,10,11. Figj2] reveals the interesting 
properties of discord. 




(46) 
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FIG. 1: The discord versus the dimensionless inverse temperature j3 for the spin chain 

with N = 21; (a) the 1st node is initially polarized (j = 1), the discord increases from the ends 
to the center of the chain; (b) the 6th node is initially polarized (j = 6), the discord is a periodic 
function of the node n. 
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FIG. 2: The discord Qj +J , 1 = 1,2, 3, 4, versus the node n for the chain of N = 21 nodes and 
different initially polarized nodes j = 1, 6, 10, 11 
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1. The discord possesses the following symmetry Q^ nn+k = Q^_ n N _ n+k , which is a con- 
sequence of the equality 

(j) = g ^N+i-n) (j) (47) 

and definition (USD. 



2. If N is odd and the middle node is initially polarized (i.e. j = pp), then the discord 
Qn n +2fe+i' k = 0,1, . . . is zero, while the discord n+2k , k = 1, . . . is zero only for 
even n, namely 

„ f 0, n = 2i, i = 1,2,... 

<n +2fc =<^ fc = l,2... (48) 

I g, n = 2z + 1, z = 0, 1, . . . 

x ((iV + 1) 2 - (4tanhf) 2 )(iV+l) 2 

2 2 ((AT + 1)2- (2tanhf) 2 ) ((JV + l) 2 - 2 (2tanh §) 2 ) 

tanhf (iV + 1 + 4tanhf) 2 (iV+ l-2tanhf) 
iVTT ° g2 (JV + 1 -4tanhf) 2 (JV + 1 + 2 tanhf) + 



y/2 tanhf (JV + 1 - 2^2 tanh § N 
log 2 



JV + 1 (JV + 1 + 2 a/2 tanh f ) ' 

Here g depends only on (3 and JV and does not depend on the particular choice of odd 
n and integer k. The reason is that 

1/2 



3K N+l 



JV+1 



JV+1 



irn 
sin — 
2 




n = 2i, i = 1,2, . . . 



, , > (49) 
2 w n = 2z + l. i = 0,1,.... 



For the case JV = 21 and /3 = 10 we obtain q m 0.0061. The direct consequence of this 
property is an equality 

+l)(2ii+2/i+l) — Q(2i2+l)(2i 2 +2l 2 +l)> ^H,h,h,h- (50) 

Eq. (1501) means that we have the system of fermions (the odd nodes) with equal discords 
between any two of them. This system is a good candidate for the quantum register. 



3. If j — 1, then the profile of the discord Q^ n+l y VZ, is bell-shaped with the maximum 
Qmax j n noc } e [ JV + 1 ~ ; ] (here [a] means the integer part of a) for discord Q^ n+l y 



If I is even, then we have # n (l)| n=[ jv+W] = gn+/(l)| n= [ jm^ ] = y iv+T cos 1(#+TJ in 
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eq.flHH). If I is odd, then we have aJDl^ N+i-^ = y jjjj cos ' q "+^ 1 ) \n-[ N+ ^- 1 ] = 
cos Ww) in eq-dlSD- In our example (/3 = 10, JV = 21) we have Q™ ax « 0.0059, 
gmo^ w 0.0058, Q^ aa; « 0.0055, QJ*™ ~ 0.0051. It is remarkable that Q nin+l) ^ for 
all possible n and I, i.e. all nodes are correlated. This is a principal advantage of the 
case j — 1 in comparison with j > 1, when the discord between some nodes is zero, 
see FigJ2} 

Finally remark, that the unitary invariant discord Q^tn+i) El as a function of the node n for 
any fixed I reproduces the shape of Q^r n +i) if J = 1 an d H> but Q G * s bigger then Q 13 . Thus, 
if j = 1, then Q^r n+ n is bell-shaped with the maximal values Q^ax equal 0.0108, 0.0106, 
0.0099, 0.0093 for I — 1, 2, 3, 4 respectively; if j = 11 and Z is odd, then the unitary invariant 
discord is the same between any neig hbors, Qjf (n+J) « 0.0028; if j = 11 and I is even, then 
the unitary invariant discord is saw-shaped with the amplitude 0.0110. 



B. Discord and concurrence in a system of c-fermions 

Now we consider the c-representation of the density matrix (see Sec JIIBI) and calculate 
the discord and the concurrence between the nth and mth c-fermions. In this case the 
/9-dependence (ft = of the discord is similar to that given in Figj2] and it will not be 
discussed here. Now the entanglement is non-zero for any iV (unlike the entanglement be- 
tween the /3-fermions) . Both the discord and the concurrence evolve in time. This evolution 
for = 10, iV = 21 and different initially polarized nodes j = 1,6,11 is shown in Figj3] 
(the discord Qn( n +i) an< ^ ^ ne concurrence Q, n+1 , between the nearest neighbors) and in 
FigJH (the discord Q^ n+2 ))- Emphasize the principal differences between the discord and 
the concurrence. 

1. The concurrence is not zero only between the nearest neighbors (except for the con- 
currence Ci 3 , which is not zero in the case of the second initially polarized node, j '• = 2; 
this concurrence is not shown in Fig 12]), while the discord between any two nodes is 
not identical to zero. For instance, the discord Q c n n+2 is shown in FigJH 

2. The concurrence is not zero only between the neighbors which are closely situated 
with respect to the initially polarized node. Thus, if j = 1, then Cf 2 , and C| 4 are 
nonzero; if j = 6, 11, then only two concurrences (namely, CJ„- ±1 ) are non-zero. 
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3. After the discord Q c nm disappears, it appears again at larger times, see Figj3j unlike 
the concurrence which does not appear again. We may refer to this phenomenon as the 
"echo" in the evolution of the discord, see FigfSh-c (where the "echo" in evolution of the 
discord Q^( n+1 ) is shown for different initially polarized nodes (j = 1, 6, 11)) and FigJH 
(where the "echo" in evolution of the discord Q c n ^ n+2 ) is shown). This phenomenon is 
not observed for the concurrence. 

4. FigsJ3h-c and |H demonstrate that the discords Q n ( n+ i) and Q( n +i){n+i+i) ar e both non- 
zero over some common time interval, which is not observed for the concurrence, see 
Figj3H-f. This is especially evident for the short time intervals, t < 10. This property 
of the discord might be used, for instance, to organize the relay of the quantum discord 
from the given pair of nodes to the neighboring pair and so on. 

C. Discord and concurrence in a system of spin- 1/2 particles 

Now we consider the basis of eigenvectors of the operators Ij Z and study the quantum 
correlations between the nth and mth spins (see Sec lII CI for the reduced density matrix 
representation). In this case the concurrence is non-zero and both the discord and the con- 
currence evolve in time. Emphasize that both the concurrence and the discord are non-zero 
only between the nearest neighbors (because the discord and the concurrence in the system 
with the diagonal density matrix are zero and p^f is non-diagonal only if m — n + 1, see 
Sec lII CI) . Moreover, since the reduced density matrix for the nearest neighbors coincides 
with that for the c- representation (see Sec lII CI) , both the discord and the concurrence co- 
incide with those obtained in Sec JIVB| i.e Q s ^ +1 = Qn, n +i an d C^n+i = ^n,n+v see FigGil 
The basic features of the quantum correlations are following. 

1. Similar to the c- representation, the concurrence is non-zero only between the neighbors 
which are closely situated to the initially polarized node. Thus, if j = 1, then C±2, C23 
and C34 are nonzero; if j = 6, 11, then only two concurrences (namely, Cjj±i) are not 
zero. 

2. Similar to Sec jIV B| there is so-called "echo" in the evolution of the discord, while this 
phenomenon is not observed for the concurrence, see FigJHJ 



16 




Dimension/ess time, D t Dimensionlesx time, D t 



FIG. 3: The time evolution of the discord Q c n n+1 (a — c) and the concurrence n+1 (d — f) for 
N = 21 and different initially polarized nodes j = 1,6, 11. The insets illustrate the "echo" in the 
evolution of quantum discord 

3. A principal disadvantage of this matrix representation in comparison with the c- 
representation is the zero discord between remote nodes. 

V. CONCLUSIONS 

For the spin-1/2 open chain of iV > 10 nodes with the XY Hamiltonian we demonstrate 
that the behavior of the quantum correlations crucially depends on the basis where these cor- 
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FIG. 4: The time evolution of the discord Q n ,n+2 for N = 21 and different initially polarized 
nodes j = 1, 6, 11 

relations are calculated. Thus, both the discord and the concurrence in the /^-representation 
do not evolve in time and concurrence is zero in this case. It is remarkable that the discord is 
nonzero for the /3-representation reflecting the quantumness of the system. It is shown that 
the system of /3-fermions with the middle initially excited node possesses the subsystem of 
nodes with all equal pairwise discords (the system of odd nodes). If the first node is initially 
excited, then all nodes are correlated, but the measure of correlations (discord) is different 
(but non-zero) for each particular pair. 

The behavior of the discord and the concurrence in another fermion representation (the 
c-representation) significantly differs from their behavior in the /3-representation. Both the 
discord and the concurrence evolve in time, but the concurrence is nonzero only between 
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the nearest neighbors (there is the single exception for Cf 3 in the chain with the second 
initially polarized node), unlike the discord. In addition, the "echo" is observed in the 
discord evolution, i.e. after the discord Q n , m (with any m > n) disappears, it appears again 
at larger times. In this representation we observe the relay of discord from the first pair of 
the neighboring nodes to the next pair and so on. 

Finally, the behavior of the concurrence between the nearest neighbors in the basis of 
eigenvalues of the operators Ij Z coincides with that for the c-representation. The discord in 
this basis coincides with the discord in the system of c-fermions only for the nearest neighbors 
and equals zero (along with the concurrence) between the remote nodes (i.e. m > n + 1). 

The proposed analysis of the quantum correlations shows that the representation of a 
quantum system in the basis of eigenvalues of the operators Ij z is less applicable in quantum 
devices in comparison with both fermion representations. Since the j3- and c-representations 
exhibit more variety of quantum correlations, it is more profitable to use the fermions (in- 
stead of the real spin- 1/2 particles) as nodes in quantum devices. This example demonstrates 
that the "virtual particles" may reveal some hidden usefull properties of a given quantum 
system. 

This work is supported by the Program of the Presidium of RAS No. 8 "Development of 
methods of obtaining chemical compounds and creation of new materials" . 
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